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LONGITUDINAL VIBRATION CHARACTERISTICS OF 
1/10-SCALE APOLLO/SATURN V REPLICA MODEL 


By Larry D. Pinson and H. Wayne Leonard 
Langley Research Center 

SUMMARY 

A 1/10-scale replica model of the Apollo/Saturn V launch vehicle has been tested 
and analyzed to determine its free-free, longitudinal vibration characteristics during 
first-stage burn. A lumped -parameter analysis was employed to assess the applicability 
of such an analysis to a complex Saturn V type of launch vehicle and to investigate the 
effects of tank orthotropy and the method of representing the mass of contained liquids 
on the calculated natural modes and frequencies of the system. An approximate theory 
which accounts for tank wall orthotropy is developed for the calculation of springs asso- 
ciated with liquid-tank systems. A simple analysis is presented which shows that the 
actual value of the mass of the contained liquid may not be the correct value for incor- 
poration into a lumped-parameter analysis. Model test results are presented for first- 
stage propellant loading conditions corresponding to several time points during first- 
stage burn. 

The experiment and analysis show generally good correlation in both frequency and 
deflection shape for structural modes; however, discrepancies are found when the system 
responds in either propellant modes or in coupled structural -propellant modes. The 
results also show that vehicle structural modes are less sensitive to shell orthotropy 
than are the modes involving significant liquid-tank interactions. Also, use of an ana- 
lytical equivalent liquid mass value for contained liquids produces better agreement with 
experiment for the liquid-tank modes whereas actual liquid mass values produced better 
agreement for structural modes. 


INTRODUCTION 

As the size and complexity of space-flight launch systems have increased, the 
acquisition of reliable vibration data for use in control-system and loads -alleviation 
studies, payload environmental parameters, and Pogo instability studies has become 
increasingly difficult and expensive. Such data are usually obtained in two ways - by 
analysis and by ground tests of a full-scale prototype. The full-scale ground testing of 
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large boosters can involve large time and money expenditures and can create technical 
and logistics problems of such magnitude as to render this approach unattractive. Fur- 
thermore, data from full-scale tests must be preceded by design and construction of the 
prototype vehicle; hence, the data are not available until late in the vehicle development 
cycle. This late availability may reduce the usefulness of the full-scale data since many 
of the design and control parameters may have already been fixed. It is desirable that 
relatively inexpensive structural dynamics data be available early in the vehicle- 
development cycle in order to evaluate structural design proposals, to identify potential 
vibration problems, and to provide some reference for the evaluation of structural anal- 
ysis techniques. 

As one possible solution to the problem, the Langley Research Center has developed 
the concept of utilizing test results from detailed subscale models as a source of experi- 
mental vibration data. (See, for example, refs. 1 to 7.) Application of the dynamic - 
model concept has been further extended to the Saturn V vehicle. A 1/10-scale replica 
model, described in reference 8, which is essentially a reduced-scale replica of the pro- 
totype, and a 1/40-scale dynamic model, described in reference 9 in which the mass and 
bending stiffness properties are scaled, have been constructed and have been subjected 
to extensive ground vibration surveys in which natural frequencies, deflection shapes, 
and damping values have been determined for various boundary restraints, simulated 
flight times, and directions of excitation. Companion analytical studies have also been 
conducted, the correlation between experimental results and analytical predictions being 
used to evaluate and improve the analyses. 

The purposes of this paper are to document the experimental free -free longitudinal 
vibration characteristics of the 1/10-scale Apollo-Saturn V model for the configuration 
and weight conditions corresponding to first -stage burn as part of the Saturn V model 
program and to present a lumped -parameter analysis which includes the effects of stiff- 
ening elements such as stringers and ring frames by basing spring -constant calculations 
on equivalent, orthotropic membrane properties of the model shell structure. Experi- 
mental response frequencies, deflection shapes, and damping values are presented and 
compared with natural frequencies and mode shapes predicted by the lumped-parameter 
analysis. 


SYMBOLS 


A area 

a radius of cylinder 
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b 

b l 

C ij 

c 

d 

E 

F 

f 

h 

K 

k n 

l 

m 


m i 

n = b/a 
n l = b l/ a 


depth of lower bulkhead dome 
depth of upper bulkhead dome 
element of orthotropic proportionality matrix 

distance from liquid surface to upper bulkhead-cylinder intersection 

deflection 

Young's modulus 

force 

frequency 

height of liquid 

spring constant 

spring constants where n = 1,2, ... 32 

shell length 

mass, also c/a 

masses where i = 1,2, ... 29 


P 


external load applied to shell 


P = 1 / a 
q= h/a 
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radius of cone at its base 


radial coordinate 
meridional coordinate 
thickness of shell 

displacement of shell in meridional direction 
displacement of liquid in longitudinal direction 

displacement of upper bulkhead -cylinder intersection in longitudinal 
direction 

displacement of liquid center of gravity in longitudinal direction 
volume of liquid 
weight of liquid 

displacement of shell normal to surface of shell 

displacement of liquid in radial direction 

longitudinal coordinate 

element of flexibility matrix 

unit weight of liquid 

deflection in x-direction 

strain 

damping ratio 


Poisson's ratio 


p 


mass density of liquid 


cp 0 cone angle 

a; natural circular frequency 

Subscripts: 
a analytical 

BH bulkhead 

c cylinder 

e experimental, also equivalent 

f fluid 

i isotropic 

s meridional 

0 orthotropic 

6 circumferential 

Matrix notations: 

[ ] denotes square matrix 

1 I denotes column matrix 

Dots over symbols denote derivatives with respect to time. 

APPARATUS AND TEST PROCEDURE 
1/10 -Scale Model 

The 1/10-scale model of the Apollo/Saturn V vehicle consists of three model 
booster stages, the model payload, and model launch escape system. The model 
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components are referred to by full-scale nomenclature: for example, S-IC, S-II, and 
S-IVB stages; instrument unit (IU); Saturn LM adapter (SLA); lunar module (LM); ser- 
vice module (SM); command module (CM); and launch escape system (LES). A complete 
description of the model, scaling concepts, solutions to typical problems encountered in 
design and fabrication, and illustrative examples of subscale components are presented 
in reference 8. Therefore, only a brief illustration of model components is included in 
this paper. The complete model is shown photographically in figure 1 and schematically 
in figure 2. Also noted on the schematic are the major structural joints and subassem- 
blies. The model is 39.6 inches (1.01 meters) in diameter at the base and is 435.5 inches 
(11.06 meters) tall. Propellants and oxidizers are shown symbolically in their respec- 
tive tanks with flight-vehicle nomenclature. Simulation of these liquids in the model is 
discussed in a subsequent paragraph. 

Some details of the model construction are shown in figures 3 to 8. Figure 3 is a 
view of the base of the model with the engines removed. Some details of interest include 
actuator support assemblies for the outboard engines, engine fairings, S-IC fuel-tank 
lower bulkhead, thrust pads, and the cross-beam support for the center engine. Also 
noted is the point at which the excitation force was applied for the tests reported herein. 
Figure 4 shows the simulated LOX suction duct and a sketch illustrating installation of 
the duct in the model. Of particular interest here are the bellows assemblies in the 
lower section of the line. Resonances associated with the deflections of the simulated 
pressure-volume control (PVC) on the bellows springs is discussed subsequently. Fig- 
ure 5 is a photograph of the complete S-II stage with major substructures identified and 
figure 6 shows a closeup view of one of the simulated S-n engine assemblies. The 
accelerometers shown are oriented in the pitch and yaw planes and were monitored 
during the tests to determine the magnitude, if any, of lateral engine responses. Fig- 
ure 7 shows the complete S-IVB stage and is presented solely for the convenience of the 
reader. The SLA and the IU are shown in figure 8 with the model LM mounted inside 
the SLA. 


Liquid Propellants 

The onboard liquids, with the exception of the LH 2 , were simulated with water. The 
actual specific gravities of these liquids are 0.80 for RP-1 and 1.14 for LOX. Proper 
control of tank water levels, as discussed in a subsequent section, then permitted the 
representation of the correct propellant weight but resulted in some deviation in mass 
distribution from that of the full-scale vehicle along the length of the model. 

The low specific gravity of LH 2 dictated an extremely lightweight simulant having a 
specific gravity of 0.07. Small, hollow, styrene plastic beads with a specific gravity of 
0.07 were procured and used as the LH 2 simulant. All tests were conducted with the tank 
ullage pressure at 10 psi (6.89 x 10^ N/m2). 
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Suspension System 


The suspension system utilized to simulate the free-free flight condition consisted 
of four, symmetrically located, vertical, steel cables extending from overhead beams to 
a cradle at the base of the model. Each cable was 403 inches (10.24 meters) long and 
made of steel strands of 7 by 19 construction. The total weight including turnbuckles 
and end fittings was approximately 9 pounds (39.6 N) per cable. Figures 9 and 10 show 
the support cradle and illustrate the method of attachment to the model. Figure 9 is an 
oblique view of the cradle and model base as seen from below. The cradle is esentially 
a square frame constructed of aluminum tubing with I-beam projections extending diago- 
nally outward from the corners of the frame. Also visible in figure 9 are the center 
engine gimbal through which the excitation force was applied, RP-1 tank fill line and 
liquid-level sight tube, S-IC thrust structure cross beams, and an inside view of one of 
the outer engine fairings. Figure 10 shows a horizontal view of the base of the suspended 
model. Visible in this figure are the adjustable clevises used to connect the cradle to 
one of the four holddown posts, the cable attachment to the cradle, the fill tube and sight 
glass for the RP-1 tank, and two of the support cables. Also visible is the single-drive 
dual-post jack used to position the model during cable hookup. The four support posts 
are used to support the model when it is not attached to the cables. The outboard engines, 
present during the tests, were not attached when the photos were made. The weight 
(18 lb (80 N)) of the cradle and shaker armature is approximately equal to the weight of 
the center engine which was removed during the tests. The length, diameter, and thus 
the cable spring constant were chosen so that the highest rigid-body natural frequency of 
the model on the cables (5.7 cps at S-IC burnout) was much lower than the lowest struc- 
tural resonant frequency (30.9 cps). 

Test Weights and S-IC Propellant Levels 

Tests were performed for S-IC propellant loadings corresponding to five time points 
during the first-stage burn. The amount of liquid in each S-IC tank was scaled on a 
mass basis from full-scale values. Table I is a listing of the nominal weight conditions 
in terms of the ratio of liquid mass in the S-IC tankage to the S-IC liquid mass at lift-off 
(in percent). The table shows the weight of the liquid in each S-IC tank, the model sta- 
tion number of the liquid surface in each S-IC tank, and the total model weight for each 
time point. The percentage values listed are nominal values. Actual time points were: 
146 seconds (burnout); 110 seconds (24 percent); maximum dynamic pressure (48 percent); 
40 seconds (73 percent); and lift-off (100 percent). 
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Instrumentation 


Instrumentation for the tests consisted of four basic groups: input force generator, 
acceleration response transducers, monitoring and recording equipment, and on-site data- 
analysis equipment. The input -force generator was a servocontrolled electromagnetic 
exciter acting through a force gage at the gimbal point of the S-IC center engine. The 
servocontrol maintained a constant, preselected, root-mean-square force level. The 
acceleration response transducers were Kistler model 303 servo accelerometers located 
on the model in such a manner as to permit determination of model deflection shapes in 
the longitudinal direction. A typical accelerometer attachment is shown in figure 11. 

The figure shows one of the accelerometers on the S-II aft skirt. This transducer was 
photographed during lateral tests; however, the attachment procedure is the same as that 
for longitudinal testing. A square solid-base plate was bonded to the surface of the model 
and four equally spaced holes were tapped in the base plate. An accelerometer mounting 
bracket was then attached to this plate, and finally the accelerometer was fitted through a 
hole in this bracket and fastened by bolting a retainer plate to the bracket. The equally 
spaced holes in the base plate permitted the orientation of the sensitive axis of the accel- 
erometer to be easily changed by 90°. To accomplish the reorientation, the entire 
accelerometer — retainer -plate — mounting -bracket assembly is removed from the base 
plate, rotated 90°, and then refastened. In addition, structural elements such as engines, 
simulated LM, and selected internal structural members, were instrumented. The loca- 
tion and orientation of the sensitive axis of each primary transducer is shown in the 
drawing in figure 12. The transducers are represented by solid diamond symbols, the 
acute angle of the diamond corresponding to the direction of the sensitive axis. Station 
numbers are given referenced to the S-IC gimbal plane (which is station 10.0 inches 
(0.254 meter)). 

Additional instrumentation consisting of strain and pressure transducers was selec- 
tively located in the S-IC structure. The locations, type, and orientation of these trans- 
ducers are given in figure 13. Each upper bulkhead was fitted with a pressure transducer 
and three pairs of strain gages. The gages were located along a meridian of the bulkhead 
and were oriented so that each pair of gages could detect horizontal, or circumferential, 
strains and vertical, or meridional, strains. Each pair of gages was assigned an identi- 
fication number and the letters H and V were used to denote horizontal or vertical 
sensitivity. The lower bulkheads were instrumented with pressure transducers and with 
six pairs of strain gages located on orthogonal meridians oriented to detect circumferen- 
tial and meridional strains. Each tank was also instrumented at selected tank wall loca- 
tions with longitudinally sensitive strain gages. Other strain gages were located at both 
ends of the LOX suction duct tunnel which contains the simulated LOX suction duct and 
on selected tension straps at the forward end of the S-IC forward skirt. A single pressure 
transducer was installed in the closure plate of the simulated PVC. 
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Data Acquisition and Reduction 

The data-acquisition system used in this investigation is shown schematically in 
rfigure 14. Data signal cables from each transducer were connected to data group 
switches in prearranged combinations so that a desired group of 12 selected channels of 
information could be monitored on either the oscilloscope, meter, or x,y plotter. The 
data were recorded on analog tape for subsequent reduction by digital computer. Pre- 
liminary qualitative analyses were also performed on site. 

All data recorded on analog tape were converted to digital format through a 
24 -point-per-cycle direct conversion. The data were then analyzed in order to determine 
the amplitude and phase angle with respect to the input force of the fundamental harmonic. 
The amplitude of the fundamental harmonic, when plotted as a function of the transducer 
station number, was then used to define the model deflection shapes. 

To obtain damping data, the force was removed from the model by suddenly cutting 
power to the shaker armature. The amplitude decay of a selected transducer was 
recorded as an oscillogram. The amplitudes of several cycles of the decaying signal 
were plotted and the slope of a straight line faired through these points was used to com- 
pute the damping value. To insure realistic data, at least 10 cycles of the decaying signal 
were used for each damping value. 


Test Procedure 

The experimental data were obtained in basically the same manner for all test con- 
ditions. The input force generator was programed to supply a constant -amplitude, sinu- 
soidal force to the model at slowly increasing frequency. Selected output signals from 
the various transducers were automatically plotted as diagrams of total amplitude as a 
function of frequency. Selected signals were simultaneously monitored as Lissajous 
patterns on an oscilloscope and peak frequencies and phase shifts were noted. Individual 
peak response frequencies were then examined by manually tuning the frequency of excita- 
tion to obtain peak response and recording the output of all transducers for subsequent 
harmonic analysis. The input force was then suddenly removed and the decaying signals 
were recorded as oscillograms from which the previously discussed damping values 
were measured. 


ANALYSIS 

The mathematical representation of the 1/10-scale Saturn V model consists of 
masses irq concentrated at discrete points along the vehicle's length connected by 
massless linear springs k n . Figure 15 shows this representation and table II gives the 
mass and spring values. This procedure has been used in other analyses, for example, 
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those of references 10 and 11. The procedure is based upon the assumption that the 
motion of a system with an infinite number of degrees of freedom can be approximated 
by analyzing an appropriate spring-mass system which has a finite number of degrees 
of freedom. The relationship between spring-mass and continuous systems is discussed 
in reference 12. The selection of the number of masses to be used in a particular analy- 
sis is somewhat arbitrary as is the manner in some cases in which these masses are 
interconnected by the springs. The approach used herein is to select an arrangement 
which results in a single type of geometry to analyze to obtain stiffness representations 
between two mass locations. This method locates mass concentration points where 
changes in shell geometry occur or where branch masses join the main structure. An 
alternate approach is to select displacements at closely spaced vehicle stations as gen- 
eralized coordinates, assume a displacement distribution between succeeding coordinate 
locations, and derive stiffness and mass matrices from the potential and kinetic energies 
which are calculated by using the assumed displacement functions. A notable feature of 
this method is that the mass matrix is not, in general, diagonal. This procedure is dis- 
cussed in reference 13. 

Some of the results of the present analysis have been presented in reference 14; 
however, the results presented herein reflect improvements in the lumped parameter 
analysis of reference 14. In the present mathematical representation of the 1/10-scale 
Saturn V model, an improved set of orthotropic spring constants were derived for the 
S-IC intertank section, the model LM representation was extended to include two degrees 
of freedom, and the elastic properties of the IU/S-rVB forward skirt joint were empiri- 
cally revised based upon static load -deflection tests of the model. These static tests 
and isolated LM tests are discussed in appendix A. 

The natural frequencies and mode shapes of the Apollo/Saturn V model were 
found by generating and solving the standard eigenvalue problem 

[K-Mu> 2 ]{y| = |0} (1) 

by a method of the Givens type (see ref. 15) where [K] is the stiffness matrix, [M] is 
the mass matrix, co is the natural circular frequency, and jy| is the mode shape. 

Three weight conditions were considered: (1) 100 percent with the liquid oxygen 

and fuel tanks of the S-IC stage filled to simulate the lift-off weight condition, (2) 50 per- 
cent, corresponding roughly to the state at maximum dynamic pressure, and (3) empty, 
corresponding approximately to S-IC burnout. 
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Mass Distribution 


The concentrated masses were found, in general, by locating one -half the structural 
.imass between the end points of a shell at each of the end-point stations. A shell is taken 
here to mean a segment of the main structure which is essentially comprised of a single 
type of geometry, such as an intertank section or an interstage adapter. If one of the 
end-point stations was the junction point of a tank bulkhead and tank wall, the entire mass 
of the bulkhead was concentrated at the junction station. Engines in the upper stages 
were considered as separate concentrated masses. The outer engines in the first stage, 
however, were grouped with the support cradle and one-half the aft skirt mass. This 
arrangement was used because the attachment of these engines to the main structure is 
virtually rigid. The simulated lunar module (LM) was considered to be a two -degree -of - 
freedom branch system. Each liquid element in the vehicle was represented by a single 
lumped mass. The plastic beads, used to simulate the liquid hydrogen in the S-n and 
S-IVB stages, were treated in the same manner as the water which was used to simulate 
the liquid oxygen in all the stages and the fuel in the S-IC stage. 

Because the liquid motion is potentially more complex than that of any other com- 
ponent of the vehicle, care must be exercised in its representation. A representation 
which has been used is that which is proposed by Wood in reference 16 and which is 
applied, with slight modifications, in references 10 and 11. The liquid is represented as 
a mass supported with certain springs and the coordinate used to describe the motion is 
the deflection of the liquid center of mass. The magnitude of the concentrated mass is 
equal to the total liquid mass. This representation was used for the liquid masses given 
in table II. A discrepancy between the analysis and experiment, which is discussed in 
the section "Results and Discussion," however, led to further consideration of this repre- 
sentation of the liquid. 

The following analysis shows that the actual liquid mass is not necessarily appro- 
priate for use with the center -of -mass deflection when this deflection is selected to be 
the generalized coordinate which characterizes the liquid motion. The kinetic energy of 
the liquid is written in terms of assumed motions. The magnitudes of these motions are 
determined by the deflection of the liquid center of mass. This deflection is therefore 
used as a generalized coordinate. The equivalent mass associated with this generalized 
coordinate is derived from the kinetic energy; thus, a mass consistent with the Rayleigh - 
Ritz procedure results. The analysis thus indicates approximately the deviation of the 
equivalent mass from the actual mass when the center -of -mass deflection is chosen as 
the generalized coordinate. 

Consider the greatly simplified propellant tank of figure 16. It is assumed that the 
bending stiffness of the cylindrical shell which contains the liquid is negligible. With 
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this assumption, the radial displacement of the shell under static conditions is 



( 2 ) 


where y is the unit weight of the liquid, E is Young's modulus, t is the shell thick- 
ness, and a and x are defined in figure 16. 

If it is assumed that a disk-shaped element of the liquid remains disk-shaped after 
the tank deformation, the longitudinal displacement of the liquid, denoted by u, is given 
by 


u = 2 



(3) 


where | is a dummy variable of integration. Expressing w in terms of uj, the 
movement of the center of mass of the liquid, gives 


u = -^u^h 2 - x 2 ) (4) 

2h 2 1V 

It is assumed that the radial displacement of the liquid, denoted by w, varies linearly 
from zero at the center of the tank to a value equal to the shell displacement at the tank 
wall. The radial displacement of the liquid is then given in terms of u^ by 


w=2r|ui 
2h 2 1 


(5) 


The kinetic energy of the liquid can be expressed in terms of Uj which is now regarded 
as a generalized coordinate. The expression for the kinetic energy T is 


T = 



dr dx 


( 6 ) 


The equivalent mass m e is found by differentiating the kinetic energy twice with respect 
to the time derivative of the generalized coordinate u^: 


m e = 


9 2 T 

^7 


From this procedure it is found that 



where m^ is the actual mass of liquid. 


(7) 

( 8 ) 
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Equation (8) shows that if the shell responds to axisymmetric dynamic loading with 
the same characteristic shape as in the case of hydrostatic loading and if the liquid motion 
is that which is assumed, the actual mass of the liquid should be increased by more than 
20 percent to obtain the equivalent mass. 

It should be pointed out that the factor derived here was obtained through the use of 
a greatly simplified model of a propellant tank. The derivation did not account for the 
presence of bulkheads and the distribution of the radial and longitudinal displacements of 
the liquid was obtained by assumption. The adequacy of these assumptions is discussed 
in the section "Propellant Tanks." 

Stiffness Representation 

The calculations made to obtain the spring constants fell into three categories: 

(1) calculations for shell components having a conical shape, of which the cylinder is a 
special case, (2) calculations for springs to represent the stiffnesses in the liquid-tank 
interactions, and (3) calculations for branch systems such as the lunar module and 
engines in the S-II and S-IVB stages. 

In order to include as accurately as practicable for a lumped-parameter analysis 
the effect of stiffening elements such as stringers and ring frames, the spring constants 
were derived by using equivalent orthotropic properties of the shells which comprise the 
total model. The major effects can be included relatively simply if the theory used is 
restricted to an axisymmetric membrane theory for orthotropic shells. This theory, 
when isotropic elastic constants are used and the shell is a cylinder, produces the con- 
ventional spring constant of AE/Z. As a consequence of the neglect of bending stresses, 
all eccentricity effects are also ignored. Reference 17 is devoted exclusively to aniso- 
tropic shell theory and it includes the orthotropic membrane theory used herein. 

Conical shell .- A typical conical shell and its associated coordinate system are 
shown in figure 17. The coordinate in the meridional direction is s and the coordinate 
in the circumferential direction (not shown) is 6 . 

The spring constant for the conical shell is defined as 


where F is the force and 6 is the distance in the x-direction through which it acts, 
if it is assumed that the base is not allowed to move longitudinally. 

The relation between the stress resultants and the strains for a membrane shell 
of revolution loaded axi symmetrically is, in matrix form, 
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( 10 ) 



C 11 C 12 


C 12 C 22 


where N s is the stress resultant in the meridional direction, is the stress resul- 

tant in the circumferential direction, e s and eg are the corresponding strains, and the 
matrix pjjj is a proportionality matrix which reflectes the orthotropic elastic proper- 
ties of the shell. 

For the force F, the stress resultants are found from equilibrium to be 


N s =- 


N 0 =O 


27t(R 0 - l cot cp Q + s cos <p 0 )sin <p Q 


(ID 

( 12 ) 


The strains are, from reference 18, 


e s = 


au 

8s 


u cos cp Q + w sin cp Q 


6e 'R. 


I cot (p 0 + S COS cp 0 


(13) 

(14) 


where u is the displacement in the meridional direction and w is the displacement 
normal to the shell. 

When equation (10) is inverted and the previously determined relations for the stress 
resultants in terms of the load (eqs. (11) and (12)), and the strain displacement relations 
(eqs. (13) and (14)), are substituted, two simultaneous equations result which may be used 
to determine displacements in terms of the load. These equations are: 


8u 

9s 


271-Cjj 1 1 - 


C 12 C 12 ' 

C 11 c 22. 


(R 0 - l cot cp Q + s cos (p 0 ^sin cp Q 


w 


= -u cot cp Q + 


C 12 F 


c 10 c, 


27TC..C,,,, 1 - -12 sin 2 <p 0 

11 22 l C n C 2 J 


(15) 


(16) 
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When equations (15) and (16) are solved simultaneously and the relation 


5 = u(0)sin <p Q - w(0)cos cp Q 


( 17 ) 


is applied, the displacement is given by 

F 


6 = 


2 7t c h C 22 f 1 -^^) sin2 < 
C 11 ^22 


(c 99 log — C 19 cos 2 <p ] 

V 22 e R o-l cot cp Q 12 °j 


(18) 


<P 0 cos cp Q 

and therefore, from equation (9), the equivalent spring for a conical shell is given by 


K = 2ttC 
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1 - 


C 12 C 12 ’ 
C 11 C 22 j 


sin ^<p Q cos cp Q 


log, 


R r 


e \R 0 - l cot cp Q J C 


— cos 2 <p 0 
22 


(19) 


The equation for the equivalent spring constant for the cylinder can be obtained as 
a special case of the cone by allowing cp Q to approach the value tt/2. When this is 
done, equation (19) becomes, in the limit, 



C 12 C 12 \ 
C 11 C 22 / 


( 20 ) 


The quantity K c , which represents the equivalent spring constant for an orthotropic 
membrane cylinder, depends primarily upon and to a much lesser extent on C 22 

and This dependence suggests that shells which do not contain liquids, such as 

intertank and interstage sections, are probably adequately represented by springs deter- 
mined by the value AE/l where A represents the total cross-sectional area of the 
shell. Although this approach would yield some savings in time, some degradation in 
accuracy could result. 

An illustrative example of the calculation of an equivalent spring for a conical shell 
is given in appendix B along with the calculation for one of the engine springs. The pur- 
pose of this appendix is to show the general degree of detail employed in the analysis of 
the Apollo/Saturn V model. 

Propellant tanks.- A typical propellant tank is shown in figure 18. It is assumed 
that the behavior of the cylindrical portion of this tank can be described by equation (10); 
that is, the membrane theory for orthotropic shells is valid. The ellipsoidal bulkhead 
shown, of which the sphere is a special case, is assumed to be isotropic and to have a 
constant thickness. The height of liquid in the tank is assumed to be greater than the 
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length of the cylinder; that is, the liquid is standing within the upper bulkhead as shown 
in the figure. 

It is further assumed that for practical tanks, the contribution of the upper bulkhead 
deformation to the total center -of -mass motion is small relative to that of the cylinder 
and lower bulkhead. In the analysis which follows, the upper bulkhead deformation is 
therefore neglected. 

It is arbitrarily assumed that a structural mass is to be located at the junction of 
the upper bulkhead and the cylinder and that motion of the liquid is characterized by the 
movement of the center of mass of the liquid, that is, the liquid is represented as a mass 
supported in some manner by springs. The base of the cylinder is assumed to be fixed 
against longitudinal movement. Consistent with these assumptions, the two coordinates 
describing the tank motion are taken to be u Q and u^, as indicated in figure 18. 

The system is assumed to be linear. Therefore, the deformation is proportional 
to the load. If the loads P and W are applied as shown in the figure, the relationship 
between the load and the deformation can be written in matrix form as 



( 21 ) 


The matrix £of] is the flexibility matrix of the system. The inversion of this matrix 
yields the stiffness matrix for the propellant -tank system. It is more convenient to pro- 
ceed with the loads P and W regarded as known quantities and obtain the flexibility 
matrix than to assume that the displacements are known and obtain the stiffness matrix 
directly. 

By assuming a hydrostatic pressure variation for the liquid and neglecting any tank 
pressurization as a static load whose effect may be added subsequently to obtain the total 
stress, the equilibrium equations for the cylinder are: 
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P 
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where T = 3q + 2n - 
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and, referring to figure 20, 
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Substituting equations (22) and (23) for the stress resultants, along with equations 
(13) and (14) for the strains into equation (10), the stress-strain law, letting cp Q = and 

inverting yields 



(24) 


Integration of the first of equations (24) and application of the boundary condition u(h) = 0 
yield the displacements in terms of the applied loads. 
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When the radial displacements of the cylinder and bulkhead are known, it is possible 
to calculate the longitudinal displacement of the center of mass of the liquid. This dis- 
placement is given, as shown in reference 18, as 


_ AM 

1 = T" 


(27) 


where AM is the first moment of the total volume change of the tank about an axis 
located in the surface of the liquid. This volume change is made up of two parts — the 
volume change which occurs because of the expansion of the cylinder wall and the volume 
change which occurs because of the expansion of the bulkhead. For small deflections the 
first moment of this volume is obtained by the relation 


AM = 2ttzl J wx dx + 222^- j^n,^) + 2qG(n,i/) + q 2 F(n,i/)] 


(28) 


where y is the unit weight of the liquid and F(n,t>), G(n,v), and H(n,i/) are the func- 
tions obtained in reference 18 in connection with equivalent bulkhead stiffness. 
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If the integration of equation (28) is performed and the result is substituted into 
equation (27), equation (21) which contains the flexibility matrix may be written explicitly, 
by using the fact that u Q = u(c), as 



(29) 


where 


Kgg — 2nEt 




9 |H(n,i/) + 2qG(n,^) + q 2 F(n,j/|] 


(30) 


The matrix relating the loads to the displacements in equation (29) is the flexibility 
matrix of the propellant -tank combination. When this matrix is inverted, the resulting 
stiffness matrix is seen to represent a system having three springs. The off-diagonal 
term represents a spring connecting the two masses while the difference between the off- 
diagonal and the diagonal terms represents a spring connecting a mass to the ground. 
Application of these facts leads to the conclusion that the set of springs represented by 
the stiffness matrix is that shown in figure 18. The value of each of these springs is 
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where 


» _ 12 q 3 ~ m 3 _ C 12 C 12 q 2 - m 2 2 77 C 22 L _ C 12 C 12 \ 

pT 2 C 11 C 22 p 2 r 2 P k BH \ c ll c 22/ 

and is given by equation (20). 

If a similar analysis is performed with the assumption that the liquid does not 
stand within the upper bulkhead, the equations for Kj, K^, and Kg are the ones 
obtained by setting m = 0 in equations (31), (32), and (33). In this case q = p. 

Although the springs of equations (31) to (33) are interrelated through the various 
geometric parameters, the orthotropic constant to which each spring is most sensitive 
can be determined by an examination of each spring. The quantity Kg, as shown by 
equation (20), depends primarily upon Cji, the longitudinal membrane stiffness. Exam- 
ination of equation (31), the expression for Kj, shows that there is a product of K c 
with C22/C11. The of K c cancels the in the denominator; thus, Ki is 

shown to be most sensitive to the value of Cg2> the circumferential membrane stiffness. 
Since Ki is directly connected to the liquid mass, as shown in figure 18, it is of pri- 
mary importance in determining the natural frequency of the liquid-tank system. This 
examination also shows that if an isotropic membrane theory is used in the liquid-tank 
problem and the circumferential and longitudinal stiffnesses are substantially different, 
significant response errors could result. Examination of equation (32), the expression 
for K2, shows that this spring is most sensitive to the value of Cj2, the constant which 
couples longitudinal stress resultants and strains to those in the circumferential direc- 
tion. This quantity therefore produces an effect in orthotropic membrane shells which 
is equivalent to that of Poisson’s ratio in isotropic membrane shells. A similar exami- 
nation of equation (33), the .expression for K3, shows that its dependence is primarily 
upon Cjj, the longitudinal membrane stiffness. 

Various special cases can be obtained from equations (31), (32), and (33). The 
springs for the simplified tank of figure 18, which was analyzed in connection with the 
equivalent mass of the liquid, are obtained by setting m = n = n^ = 0, p = q, and Kgjj 
approaching °°, and by assuming isotropy. The result is 
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If the structural mass is taken to be zero, the system reduces to a single -degree - 
of -freedom system with an equivalent spring 


K e = 


37rEt 

2q 


(37) 


With this result the natural frequency of the system can be calculated by using the equiva- 
lent mass of equation (8) and the result can be compared with available data on this con- 
figuration. References 19 and 20 contain such data based upon more rigorous, but more 
complex, analyses. Figure 19 shows a comparison of the frequencies from each of these 
analyses for various tank height-radius ratios. The results show a maximum deviation 
of about 6 percent for the range of heights investigated, the greatest difference occurring 
for the most shallow tank. It is felt that these comparisons tend to justify as being rea- 
sonable the various assumptions made to obtain the greatly simplified analysis from 
which the equivalent mass of equation (8) was derived. 


RESULTS AND DISCUSSION 


Experimentally measured and analytically determined longitudinal, free -free 
vibratory characteristics of the 1/10-scale Apollo/Saturn V model are presented in fig- 
ures 20 to 51 and in table III which lists both analytical and experimental resonant fre- 
quencies and experimental damping values, where obtainable, for various S-IC liquid 
loading conditions. A summary plot of the experimental resonant frequencies as a func- 
tion of S-IC loading condition is presented as figure 20. The closed symbols of figure 20 
represent data from table m whereas the open diamond symbols are values obtained from 
supplementary, intermediate tests conducted during the course of the program, but other- 
wise not documented in the present paper. The purposes of this figure are to give the 
reader a visual summary of the model behavior and to illustrate the coupling between the 
S-IC LOX mode and structural and branch mass resonances. The faired curves repre- 
sent the authors' interpretation of the model behavior based upon nonrecorded interme- 
diate tests and recorded intermediate test points such as those shown. 

In general, the experimental data are presented in the same format for each pro- 
pellant loading condition. Diagrams of acceleration as a function of frequency are shown 
wherein the magnitude of the total response of selected accelerometers is displayed as a 
function of the forcing frequency for a constant-magnitude sinusoidal -input force. The 
trace associated with a given transducer has its zero at or near the location of the trans- 
ducer referenced to the schematic of the model as shown on the left-hand side of each of 
the diagrams. Plus and minus signs indicate, respectively, in-phase or out-of -phase 
response with respect to the input force. Although all traces on a given diagram are 
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plotted with approximately equal sensitivities, the sensitivities are not necessarily the 
_same on different diagrams. 

After the acceleration-frequency diagram, the experimental deflection patterns for 
constant -frequency excitation, and where applicable, the mode shapes calculated by uti- 
lizing the orthotropic analysis are shown. The liquid mass in these calculations was 
taken to be the actual mass of the contained liquid. After the data presentation, a discus- 
sion of analytical results using both isotropic and orthotropic membrane theory for stiff- 
ness and using the actual and equivalent mass concepts for the liquid is presented. Both 
the experimental and analytical results are normalized to unity at the forward end, or tip, 
of the model wherever possible; however, in some instances, low experimental tip ampli- 
tudes necessitated the selection of alternate normalization points. The normalized 
amplitudes are plotted as a function of the nondimensional longitudinal coordinate x/L, 
where x is the model coordinate and L is the total length of the model. On each plot, 
the experimental resonant frequency, the appropriate analytical natural frequency, and 
the measured free-decay damping value are tabulated. 

Lift-Off Configuration 

Data from the lift-off configuration (100-percent S-IC propellant loading) are shown 
in figures 21 to 30. Figures 21 and 22 show two acceleration -frequency diagrams for 
this particular configuration. The trace-deflection sensitivities are the same for the two 
plots; however, the data of figure 21 were obtained with a 6 -pound (27 N) nominal force 
input and the data of figure 22 were obtained with a nominal 12 -pound (53 N) input force. 
The data show that doubling the input force approximately doubled the response amplitude 
and thus indicated that the structure is nearly linear for this range of input force. The 
acceleration-frequency diagrams also reveal that the response of the model is divided 
into two distinct frequency domains. At lower frequencies (f e < 80), the structure above 
the S-II stage is very active whereas the remainder of the model is relatively quiescent. 
Conversely, in the higher frequency domain (f e > 80), the S-IC and S-EI structures respond 
to the input force whereas the S-IVB and payload structures experience minimal acceler- 
ations. This divided response trend was observed for all simulated flight times. 

The behavior of the trace associated with the thrust structure can best be explained 
by considering the relative locations of the input force and the accelerometer. The 
accelerometer was located on the top of the thrust structure cross beam and immediately 
above the shaker attachment point. If the cross beam is considered to be a system with 
a high stiffness-mass ratio, the accelerometer senses the reaction of a force applied to a 
mass through a spring. The lowest, experimentally determined, resonant frequency of 
the thrust structure cross beam is approximately 260 cps. Below this resonant frequency, 
the acceleration response is proportional to the square of the forcing frequency. When 
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this response is plotted as a function of the forcing frequency, as is done in the 
acceleration -frequency diagrams, a nearly parabolic curve results with superimposed 
peaks caused by system resonant responses. 

Experimental resonant deflection shapes and calculated mode shapes for the lift-off 
configuration are shown in figures 26 to 30. All data are normalized to unity at the for- 
ward end or tip of the model with the exception of figure 30 wherein low tip amplitude 
necessitated normalization at the aft end or base of the model. The analytical curves 
shown are faired through discrete calculated deflections. 

The experimental data in figure 26 show that the entire instrumented portion of the 
model structure is moving in phase with the input force. Maintenance of momentum bal- 
ance in the system precludes the existence of such a condition for the total system since 
the excitation frequency (30.9 cps) is well above the rigid-body frequency (3.5 cps) of the 
model on its support cables. With no analytical data available, it would be reasonable to 
assume that momentum balance is achieved through out-of -phase motion of some large 
noninstrumented mass. The analysis verifies this assumption and shows the out-of- 
phase mass to be the liquid in the S-IC LOX tanks. Further experimental confirmation 
of the nature of this resonance may be found from examination of the strain and pressure 
data shown in figures 23 to 25. The locations of the transducers are given in figure 13. 

In figure 23, the pressure gage P-3 located in the LOX tank lower bulkhead shows a small 
pressure response peak and phase change at about 30 cps. Likewise, in figure 24 both 
horizontal and vertical strain peaks in the LOX tank lower bulkhead are shown by the 
response of gages 22-H and 22-V. The data, however, are insufficient to determine the 
magnitude of the LOX mass motion. The surface strain data in figure 25 show no detect- 
able longitudinal strain patterns in the S-IC LOX tank at 30 cps. Further, the orientation 
of the surface gages was such that circumferential strains in the tank walls could not be 
measured. In spite of these limitations, it is felt that the experimental data along with 
the analytical results are sufficient to identify qualitatively the LOX mass resonance. 

The response of strain gages S-14 and S-16 in figure 25 also served to identify the 
LOX suction duct tunnel resonance at about 97 cps whereas the output of the pressure 
gage P-5 shown in figure 24 served to identify the resonance of the simulated PVC on its 
lower bellows spring. 

Rather poor correlation between the results of the analysis and those obtained 
experimentally is evident in both frequency (f Q = 38.0 cps; f e = 30.9 cps) and in deflec- 
tion shape. One potential reason for the poor frequency correlation is the previously 
discussed representation of the contained liquid. Effects of the application of the equiva- 
lent mass concept as derived in equation (8) are discussed in the section "Equivalent 
Mass and Orthotropic Stiffness Effects." 
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Both analysis and experiment in figure 27 exhibit the classical first bar-type longi- 
tudinal deflection shape. Frequency correlation for this mode is excellent (f Q = 40.3 cps; 
f e = 40.2 cps) and the deflection shape correlation is reasonably good with the exception 
of the simulated LM for which the experimental value is low by about 25 percent. The 
analysis also shows that the S-IC LOX mass is still moving out of phase with the tank 
structure; thus, considerable liquid-tank interaction in this mode is indicated. 

These facts illustrate the difficulties associated with the lumped -parameter model 
in treating the liquid-tank interaction. Also, they show the need for an improved model. 

A continuous analysis which was obtained by NASA contract and which is applied in ref- 
erence 14 predicts the LOX mode satisfactorily. In this analysis, several assumed 
shapes are used in the liquid-tank areas yielding more than one degree of freedom. Thus, 
work to improve lumped-parameter representations should be along the lines of including 
more degrees of freedom for the liquid. This approach is also suggested by noting that 
the equivalent mass reflects variations in the tank shape and liquid behavior which may 
occur with changes in the excitation frequency of the system; thus, a single-degree-of- 
freedom representation of the liquid will be inadequate for more than one system mode. 

The test results indicate an unusually abrupt change in the deflection shape in the 
IU-SLA joint (s/L = 0.76). This effect is present in the results for all the weight condi- 
tions. To assess the contribution of this joint to the model deformation, static tests 
were performed. The analysis was revised accordingly since the effect could not be 
derived from model drawings. These static tests are discussed in appendix A. 

The resonance depicted in figure 29 is primarily due to motion of the simulated 
LM mass. A relatively high force level was required to produce the low amplitude 
shown. For example, the tip acceleration amplitude measured for the 11-pound (49 N) 
input in figure 28 was 0.071g whereas the tip acceleration amplitude for the 15-pound 
(67 N) force in figure 29 was only 0.013g. The character of the simulated LM response 
as determined by separately testing the simulated LM is discussed in appendix A. 

The poor correlation between analytical and experimental deflection shapes in the 
S-IC tank areas in figure 29 may originate from several sources. One possible reason 
for the measured amplitudes may be harmonic excitation of the transducers on their 
mounting pads. By considering the tank wall as a rotational spring and the relatively 
massive transducer attached to the wall, cantilever resonances of the transducers on 
their mounting springs were found by plucking each transducer and monitoring the decay 
frequency. Resonances occur at approximately 200 cps for the S-IC transducers in the 
lift-off condition. All the surface transducers on the S-IC stage in the lift-off condition 
exhibited 200 cps responses with amplitudes approximately one-fourth that of the funda- 
mental at 71.3 cps. This explanation is not entirely satisfactory, however, since the 
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surface strain data of figure 25 show some evidence of longitudinal strain in the S-IC 
fuel tank walls at about 70 cps and figures 23 and 24 give indications of strain and pres- 
sure responses in the S-IC LOX tank lower bulkhead. It must, therefore, be assumed 
that the experimental deflection shape in figure 29 is at least partially associated with 
tank response patterns. The analysis does not consider higher order axisymmetric tank 
modes or nonaxisymmetric modes. The instrumentation, also, does not yield sufficient 
data to describe the response properly. 

The resonance shown in figure 30 was not predicted by the analysis; therefore, only 
experimental data are presented. In addition, there are no clear-cut indications as to the 
source or sources of the resonant amplitudes. Surface strain gages (fig. 25) indicate 
what appears to be a second LOX suction duct tunnel resonance and longitudinal strains 
in the S-IC fuel tank at 136 cps. Bulkhead strain gages (fig. 24) show both radial and 
circumferential strains in both the LOX and fuel tank lower bulkheads. S-IC pressure 
responses (fig. 23) show a large PVC pressure response at 136 cps as well as lesser 
pressure responses in both lower bulkheads. In addition, all S-IC engines and the S-IE 
thrust structure exhibit large longitudinal response peaks at this frequency. 

7 5 -Percent -Full Condition 

Experimental data from the 1/10-scale Apollo/Saturn V model with the S-IC stage 
loaded to 75 percent of its lift-off propellant weight are presented in figures 31 to 33. 

In this configuration, the lowest resonance encountered has the classical bar-type, first 
longitudinal deflection shape. The LOX mode is shifted upward in frequency and appears 
to couple with the second structural mode. The resonance of the simulated LM mass was 
barely detectable at this weight condition with nominal -input force levels. 

50-Percent-Full Condition 

Analytical and experimental data from the configuration in which the propellant 
loading in the S-IC stage is half that at lift-off are shown in figures 34 to 39. As in pre- 
vious cases, the data consist of plots of experimental acceleration against frequency and 
plots of combined analytical and experimental resonant deflection shapes. The agreement 
between experimental and analytical data for the first two longitudinal resonances of this 
condition is generally good in both frequency and deflection shape. 

It should be recognized that the experimental data presented in this paper do not 
necessarily represent a complete set of vehicle modes. Often the natural frequencies of 
a complex system are not separable by the techniques and instrumentation employed in 
this investigation since they can exist very close together. Such is probably the case in 
figure 37 wherein no experimental resonance which corresponds to the analytical reso- 
nance was found. The frequency associated with the mode presented in figure 38 
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corresponds roughly to the uncoupled frequency of the S-IC LOX mass on its supporting 
springs. However, as seen from the summary plot in figure 20, the experimental LOX 
mode is found in proximity to the LM resonance and is not separated experimentally 
from it. The coupled resonant deflection shape is presented as figure 38. This coupling 
effect as well as previously discussed bulkhead and nonaxisymmetric tank motions may 
account for the deflection shape disagreement in the region 0 S x/L ^ 0.4. The experi- 
mental resonant deflection shape of figure 39 is shown in the same context as discussed 
for figure 30. 


2 5 -Percent- Full Condition 

Experimental data from the 25-percent S-IC propellant loading condition are pre- 
sented in figures 40 to 45. As was the case with the 75-percent data, no analytical data 
were calculated. The experimental data follow the trends of previously presented 
responses with the exception that there appears to be no simulated LM response in the 
vicinity of 70 cps. Neither the acceleration -frequency diagram nor the deflection shape 
plots show evidence of any LM response at the anticipated frequency. 

Burnout Condition 

Experimentally determined response characteristics of the 1/10-scale Saturn V 
model at the S-IC burnout condition are correlated with the analytically derived mode 
shapes and frequencies in figures 46 to 51. The lowest mode (fig. 47) agrees very well 
in both frequency and deflection shape with the experimental data. Evidence of the insep- 
arability of proximate modes by conventional experimental techniques is shown in fig- 
ures 48 and 49. Analytically, two natural modes were found at 70.8 cps and 74.7 cps. 

A single resonance at 69.6 cps was observed experimentally. Examination of the exper- 
imental deflection shape reveals both the LM resonance, which is predicted by the analy- 
sis as the dominant mass in motion at 70.8 cps, and significant motion in the region 
0 § x/L S 0.5, which is characteristic of the analytical mode shape at 74.7 cps. From 
these observations, it would appear that the experimental shape represents a coupling of 
the two analytically predicted mode shapes. 

Frequency agreement between analysis and experiment for the third resonance 
(fig. 50) is only fair and there is some discrepancy in the magnitude of the deflection 
shape in the region 0 < x/L < 0.4. The experimental resonant deflection shape at 
135.1 cps is presented, previously noted comments for figure 30 being applicable. 

Equivalent Mass and Orthotropic Stiffness Effects 

The present paper characterizes as an improvement the reflection of the ortho - 
tropic properties of the structure in the calculation of system spring constants, and 
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equation (8) in the analysis indicates a potential improvement in representation of liquid- 
tank interactions by the inclusion of a liquid equivalent mass. In an effort to define qual- 
itatively the effects of both orthotropy and liquid mass value on the natural frequencies 
and mode shapes of launch vehicle structures, natural frequencies and mode shapes of 
the 1/10-scale Apollo/Saturn V model were calculated for three S-IC propellant loading 
conditions (100-, 50-, and 0-percent full) using both isotropic and orthotropic spring con- 
stants and both the actual liquid mass and the liquid equivalent mass as defined in equa- 
tion (8). For all modes except the lowest two modes in the lift-off weight condition, 
results from the orthotropic and isotropic analyses were nearly identical. Since these 
two modes exhibit the most pronounced liquid-tank interactions, it would appear that the 
orthotropy of the structure affects primarily modes of this type. This observation tends 
to support the rationale previously developed in the discussion of propellant tank stiffness. 

Four calculated mode shapes and frequencies for the lowest, or LOX, mode of the 
1/10-scale Apollo /Saturn V model are shown in figure 52. Figures 52(a) and 52(b) are 
results of the isotropic analysis using the actual liquid mass and equivalent liquid mass, 
respectively. Figures 52(c) and 52(d) are results of the orthotropic analysis with the two 
liquid mass representations. A comparison of the four figures with each other and with 
the experimental data in figure 26 shows that, as would be expected, the higher equivalent 
mass tends to lower the resonant frequency. The mode shapes associated with the equiv- 
alent mass are also distorted in the region 0 ^ x/L S 0.25. In none of the cases is the 
mode shape correlation between analysis and experiment considered good; however, the 
orthotropic theory with actual liquid mass yields reasonably good agreement. From a 
frequency standpoint, all the calculated frequencies are high, ranging from 39.3 cps to 
33.8 cps as compared with the experimental frequency of 30.9 cps. In both analyses, 
however, use of the equivalent mass yields much better frequency correlation, the ortho - 
tropic analysis yielding the closer value. 

The calculated frequencies and mode shapes for the second mode, corresponding to 
the first structural mode of the model, are shown in figure 53. As in the previous figure, 
the plots shown result from application of the orthotropic and isotropic analyses with 
actual and equivalent liquid mass representations. In this mode, the variations in both 
mode shape and frequency are less pronounced than in the LOX mode. A comparison with 
the experimental data in figure 27, however, shows that the orthotropic analysis with 
actual liquid mass yields the better mode shape correlation and nearly identical frequency 
(f e = 40.2 cps; f 0 = 40.3 cps). The degree of liquid-tank interaction in this mode as com- 
pared with the LOX mode is small. Therefore, the sensitivity of both mode shape and 
frequency to orthotropy in the tank structure is reduced. In the higher structural modes 
for this weight condition and in all structural modes for the other weight conditions ana- 
lyzed, the sensitivity to orthotropy becomes negligible. 
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CONCLUSIONS 


A 1/10 -scale replica model of the Apollo/Saturn V launch vehicle has been tested 
and analyzed to determine its free-free longitudinal vibration characteristics during 
first-stage burn, A lumped -parameter analysis was employed to assess the applicability 
of such an analysis to a complex Saturn V type of launch vehicle and to investigate the 
effects of tank orthotropy and the method of representing the mass of contained liquids 
on the calculated natural modes and frequencies of the system. In addition, an approxi- 
mate theory which accounts for tank -wall orthotropy is presented for the calculation of 
springs associated with liquid-tank systems. Model test results are presented and uti- 
lized in performing comparative evaluation of the various analytical predictions. Salient 
results of the combined analytical and experimental program may be summarized as 
follows: 

1. Launch vehicle structures in which relatively large liquid masses are contained 
in thin walled tanks exhibit three principal types of longitudinal resonant response; one 

in which the motions are primarily those of the liquid masses deforming the elastic tanks; 
one in which bar -type structural deformations are involved; and the third in which elas- 
tically mounted components respond. The model tests show that system modes involving 
significant coupling exist when the frequencies of the propellant modes occur near struc- 
tural or branch mass frequencies. 

2. By using model tests as the standard for comparison, the analysis in which the 
stiffness coefficients reflect orthotropic properties of the structure yielded significantly 
better results for the modes involving liquid-tank interactions than the isotropic analysis 
did. Comparison of the two analyses for the structural modes showed little difference in 
either the mode shapes or natural frequencies. These findings indicate that the liquid- 
tank responses are much more sensitive to orthotropy than the structural responses are. 
This result would also indicate that even for liquid-tank modes, isotropic representation 
of non -liquid-carrying shells would result in little degradation in accuracy of calculated 
modes or frequencies. 

3. Calculations performed by using both orthotropic and isotropic analyses and 
with both actual and equivalent liquid mass values show that use of an equivalent liquid 
mass produces better frequency agreement with experiment for the liquid-tank modes 
whereas the analyses utilizing actual liquid mass values were more accurate in the pre- 
diction of structural mode shapes and frequencies. The data indicate that a single- 
degree -of -freedom representation of the liquid is inadequate. 

4. A comparison of the tank -bulging-mode frequency as calculated by refined 
methods and that calculated by the approximate method developed in the paper showed 
almost identical results for a wide range of ratios of liquid depth to tank radius. 
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5. Static and dynamic tests showed an abrupt change in model longitudinal deflection 
shapes at a joint in which thin-gage materials were used. This change necessitated 
revision of the analytical stiffness for this area on the basis of static tests since this 
effect could not be derived from model drawings. It is concluded therefore that stiff- 
nesses in areas containing joints constructed of thin-gage materials should be verified 
by experiment whenever possible. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 21, 1969, 

124-08-05-18-23. 
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APPENDIX A 


SUPPLEMENTARY TESTS 

Longitudinal Static Load -Deflection Measurement 

Examination of the dynamic -response data indicated that the joint between the IU 
and the SLA (x/L = 0.76) was contributing substantially to the deformation of the model. 

An examination of the joint revealed that the mating surfaces were not planar and that vis- 
ible gaps existed on the circumference between fasteners. These gaps are attributed to 
limitations in the technology of the fabrication of thin -gage metals since the flanges of the 
ring frames which mate to form this joint are each only 0.009 inch (2.29 x 10“4 meter) 
thick. In order to determine the magnitude of this joint effect, the model was statically 
loaded and longitudinal displacements were measured along the external surface. Results 
of one such load -deflection test are shown in figure 54. The figure shows experimental 
and calculated longitudinal deflections resulting from a 400-lb (1780 -N) load applied along 
the model center line through a plate at the top of the simulated service module. Since 
there was no way to predict this joint effect from the model drawings, the spring kg of 
the analysis was empirically revised to obtain an equivalent stiffness between stations 
310 inches (7.875 meters) and 326 inches (8.277 meters). The figure shows calculated 
deflection curves from both the original and revised analyses. The representation 
employing the revised value of kg is used in all the analytical results shown in this 
report. 


Vibration Tests of Simulated LM 

Longitudinal vibration tests of the 1/10-scale Apollo/Saturn V model revealed a 
major resonance which appeared to consist primarily of motion of the simulated LM. 
Calculated spring constants for the LM body on the attachment truss members yielded 
predicted responses of much higher frequency than those observed. Consequently, the 
simulated LM was removed from the model and tested to determine the nature of the 
observed resonance. A simple test setup was employed as shown pictorially in figure 55. 
Each support truss apex was rigidly attached to an aluminum post which was, in turn, 
bolted to deep web steel I-beams set in a concrete floor. A small shaker was attached at 
the center of the aft end of the model LM. Responses were measured at selected points 
on the model and support structures. Some typical results are shown in figure 56. 
Acceleration-frequency plots from six of the transducers are given. The data indicate a 
resonance of the model LM at 67.9 cps. Examination of the vibrating model showed large 
deflections at the model center line; smaller in-phase deflections near the edge of the 
model body; no transverse motions; and no response of the support structure. It was, 
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therefore, concluded that the resonance was that of the internal ballast masses deflecting 
with the elastic model end plates. A new analysis was performed in which the LM model 
was treated as a two -degree -of -freedom system. One mass represented the ballast 
weights and end plates and the other mass represented the cylinder and attachment 
trusses. This representation predicted a LM model resonance at 67.7 cps. Because of 
the good agreement between the calculated and experimental resonant frequencies, the 
two-degree-of -freedom representation was incorporated into the mathematical model of 
the system. 

The LM model was also statically tested and the results verified the initially calcu- 
lated spring constants for the attachment truss members. 
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TYPICAL STIFFNESS CALCULATIONS 

Some typical calculations are presented to illustrate the methods used in obtaining 
the orthotropic elastic membrane constants for the shells of the Apollo/Saturn V model 
which were used in the calculation of spring constants and to show the general degree of 
complexity of the structural theory used to calculate the spring constants in the spring- 
mass model. Accordingly, the numerical calculations associated with the S-IC liquid 
oxygen tank and those for a branch mass, the S-II center engine, are presented. 

S-IC Liquid Oxygen Tank 

The S-IC liquid oxygen tank is comprised of an aluminum skin which has constant 
thickness in the circumferential direction and which tapers linearly in the longitudinal 
direction, integrally milled, longitudinal, tee-section stiffeners, and baffle assemblies 
which are welded to the flanges of the longitudinal stiffeners. All stiffening elements are 
also aluminum. Typical repeating elements are shown in figure 57. 

It was assumed that only the skin and stringers are effective in resisting the longi- 
tudinal loads and that the skin is the only element which couples longitudinal load to cir- 
cumferential deformation and circumferential load to longitudinal deformation. Since the 
skin of this tank varies in thickness from 0.019 inch (4.8 x 10 “4 meter) at the forward 
end to 0.025 inch (6.3 x 10 “4 meter) at the aft end, it was assumed that the average value 
of 0.022 inch (5.5 x 10"^ meter) could be used. Poisson's ratio is assumed to be 0.3. 
Young’s modulus is taken to be 10^ lb/in^ (6.89 x 10*0 N/m2). 

Since a portion of the baffle web is corrugated, it cannot be fully effective in 
resisting the circumferential loads. Consequently, in the calculations which follow, the 
effective area of the corrugated portion of the baffle web was taken to be 50 percent of 
the actual area. The actual area of this component constitutes about 10 percent of the 
total area of the baffle assembly and skin. Therefore, even if the 50-percent factor is in 
error by as much as 30 percent, the membrane constant would contain only a 3-percent 
error as a result of this assumption. 

For the longitudinal direction in which the total number of stringers is 168 the cal- 
culation for yields the following values: 

Stringer flange area, 168 x 0.13 x 0.015 0.32760 in2 (211.3 x 10“4 m^) 

Stringer web area, 168 x 0.163 x 0.010 0.27384 in^ (176.6 X 10~4 m2) 

Total stringer area 0.60144 in2 (387.9 x 10 - 4 m2) 

Circumference, 2 tt x 19.8 124.407 in. (3.16 m) 
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Average area of stringers per inch of circumference, 0.00483 in^ (3.1 x 10“4 m) 

c _ 0.022 x 10? + 0.00483 x 10 7 = 2.901 X 10 5 lb/in. (5.080 X 10 7 N/m) 

11 1 - (0.3)2 

For the circumferential direction in which the total number of baffle assemblies 
is 13 the calculation for C 22 yields the following values: 

Outer ring flange area, 13 x 0.015 x 0.15 0.02925 in2 

Outer ring web area, 13 x 0.145 x 0.012 0.02262 in^ 

Baffle web area (corrugated), 13 x 2.460 x 0.004 X 0.5 . . . 0.06396 in2 
Baffle web area (noncorrugated), 13 x 0.305 x 0.004 .... 0.01586 in 2 

Inner ring flange area, 13 X 0.203 x 0.011 0.02903 in 2 

Inner ring web area, 13 x 0.114 x 0.011 0.01630 in 2 

Total circumferential stiffener area 0.17702 in 2 

Shell length, 48.898 inches (1.24 meters) 

Average circumferential stiffener area per inch of length, 0.00362 in 2 (0.02 x 10"4 m 2 ) 

c _ M. 22 x 10 7 + 0.00362 X 10 7 = 2.780 X 10 5 lb/in. (4.869 X 10 7 N/m) 

22 1 - (0.3) 2 

The coupling coefficient is calculated by considering only the skin as a coupling 
element. This calculation yields: 

c _ 0.3 x 0.022 x 10 7 _ 0>7 253 x i 0 5 lb / in> ( 1>2 70 x 10 7 N/m) 

1 - (0.3) 2 


(0.19 x 10- 4 m2) 
(0.15 x 10-4 m 2) 
(0.41 x 10-4 m 2 ) 
(0.10 x 10-4 m 2) 
(0.19 x 10-4 m 2 ) 
(0.11 x 10-4 m 2 ) 
(1.14 x 10“4 m2) 


S-n Center Engine Support Spring 

The S-II center engine is supported by a cross beam of I-configuration. For pur- 
poses of analysis, the cross beam was regarded as four identical beams, each of which 
is assumed to have zero slope at the vehicle center line and to be pinned at the juncture 
of the cross beam and thrust cone. The longitudinal spring constant for the S-II cross 
beam is taken to be four times the reciprocal of the relative end deflection when a unit 
shearing force is applied to one beam. Figure 58 shows the configuration of the beam 
and the area moment of inertia of the cross section as a function of the distance from the 
center line. Shear deformation is neglected. The beam is divided into sections Axj 
1 inch (0.0254 m) long and the moment of inertia is assumed to be a constant equal to the 
average value in each of the sections. 
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The calculation of the spring constant for the cross beam follows: 


Moment, 
l(z - Xi) 

Distance to 
center line of 
section, Xj 

Average moment 
of inertia, 1^ 

lb -in. 

N-m 

inch 

meter 

inch^ 

meter^ 

8.67 

0.22 

8.67 

0.22 

0.0896 

3.73 X 10-8 

7.67 

.19 

7.67 

.19 

.2238 

9.31 

6.67 

.17 

6.67 

.17 

.1681 

6.99 

5.67 

.14 

5.67 

.14 

.1325 

5.51 

4.67 

.12 

4.67 

.12 

.1035 

4.31 

3.67 

.09 

3.67 

.09 

.0790 

3.29 

2.67 

.07 

2.67 

.07 

.0585 

2.43 

1.67 

.04 

1.67 

.04 

.0403 

1.68 

.67 

.02 

.67 

.02 

.0424 

1.76 

.085 

.002 

.085 

.002 

.0215 

.89 


10 , 

= ^ i ^ i 1 = 2192 x IQ” 7 in./lb (12 494 X 10“ 10 m/N) 


i=l 


where d is the deflection due to a unit shear on one beam. The spring constant K for 
all four cross beams is 


K = - - = 0.1824 x 10 5 lb/in. 

2192 x '10- 7 


(0.3194 x 10 7 N/m) 
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TABLE I.- TEST WEIGHTS AND S-IC PROPELLANT LOADINGS 


Nominal 

weight 

condition 

S-IC LOX 

S-IC fuel 

Total model 
weight 

(a) 

Weight 

Liquid surface 
station 

Weight 



Liquid surface 
station 


pounds 

1 . .. 

newtons 

inches 

meters 

pounds 

newtons 

inches 

meters 

pounds 

newtons 

Burnout 

4 

18 

78.3 

1.99 

12 

54 

24.5 

0.62 

1680 

7 467 

25 percent 

784 

3 488 

99.5 

2.53 

284 

1262 

33.5 

.85 

2732 

12 143 

50 percent 

1509 

6 713 

115.8 

2.94 

593 

2638 

40.5 

1.03 

3766 

16 739 

75 percent 

2274 

10 116 

133.0 

3.38 

922 

4102 

47.9 

1.22 

4860 

21 602 

100 percent 

2985 

j 

' 13 278 

151.5 

‘ 3.85 

1282 

| 5704 

56.0 

1.42 

; 5931 

26 362 


a Does not include 15.25 pounds (67.8 N) support cradle or 36 pounds (160 N) of support cables. 


TABLE n.- STIFFNESS AND MASS DATA FOR LUMPED -PARAMETER ANALYSIS 


Branch mass 

Mass and 
spring constant 
numbers 

Station 

Mass 

Spring constant 

name 

inches 

meters 

lb-sec 2 /m* 

kg 

lb/in. 

N/m 


1 

420.0 

10.688 

0.011296 

1.978 

0.1414 x 10 6 

0.2477 x 10 8 


2 

384.0 

9.754 

.019906 

3.486 

1.0572 

1.8514 


3 

377.9 

9.599 

.072545 

12.71 

.0684 

.1198 


4 

359.6 

9.130 

.066717 

11.68 

.0755 

.1322 


5 

332.7 

8.450 

.003857 

.6754 

2.0850 

.0226 

LM 

6 

332.7 

8.450 

.017090 

2.992 

.0097 

.01698 

LM 

7 

332.7 

8.450 

.05132 

8.987 

.4120 

.7215 


8 

325.9 

8.277 

.010569 

1.851 

.1220 

.2137 


9 

310.1 

7.875 

.021213 

3.715 

.3802 

.6658 


10 

283.2 

7.193 

.024752 

4.335 

.0647 

.1133 

S-IVB engine 

11 

264.6 

6.721 

.008390 

1.469 

.0969 

.1697 

S-IVB LOX 

12 

283.2 

7.193 

.45336 

79.40 

.1876 

.3285 

S-IVB LH 2 

13 

301.5 

6.758 

.104904 

18.37 

.5198 

.9103 


14 

274.7 

6.976 

.02388 

4.201 

1.4485 

2.5367 


15 

251.9 

6.398 

.024788 

4.341 

.5487 

.9609 


16 

238.7 

6.063 

.058743 

10.29 

1.2956 

2.2689 


17 

184.8 

4.694 

.074142 

12.98 

.3455 

.6051 


18 

178.7 

4.539 

.092255 

16.16 

.0641 

.1123 

S-XI outer engine 

19 

166.4 

4.227 

.033947 

5.945 

.1185 

.2075 

S-H center engine 

20 

166.4 

4.227 

.011089 

1.942 

.3916 

.6858 

S-H LOX 

21 

184.8 

4.694 

1.9933 

349.1 

5.0894 

8.9129 

s-n lh 2 

22 

211.8 

5.379 

.4059 

71.08 

1.2566 

2.2006 


23 

140.1 

3.559 

.136495 

23.90 

.0182 

.0319 


24 

91.2 

2.317 

.111689 

19.56 

.7002 

1.2262 


25 

60.2 

1.529 

.092247 

16.16 

a .6162 

1.0791 


26 

36.5 

.927 

.178389 

31.24 

a o 1008 

.1765 


27 

11.2 

.285 

.40678 

71.24 

a .2752 

.4819 

S-IC fuel 

28 

41.3 

1.050 

3.3294 

a 583.1 

.8584 

1.5033 

S-IC LOX 

29 

115.7 

2.938 

7.6627 

a 1342. 

1.2056 

a 2.1113 


30 





.0778 

a 1.362 


31 





a .4902 

a .8585 


32 


1 



3.1748 

5.5599 


a The values shown apply to the 100-percent condition. In the 50-percent condition: m 28 = 1.6647 (291.6); 
m 29 = 3.8313 (671.0); k 25 = 0.6351 x 106 (1.1122 x 10 8 ); k 26 = 0.0604 x 10 6 (0.1058 x 10 8 ); k 27 = 0.6000 x 10 6 
(1.0508 X 10 8 ); k 29 = 1.2622 x 10 6 (2.2105 x 10 8 ); k 30 = 0.0108 X 10 6 (0.0189 x 10 8 ); and k 31 = 0.5568 x 10 6 
(0.9751 X 10 8 ). In the burnout condition: m 28 = m 29 = 0; k 25 = 0.6900 x 10 6 (1.2084 x 10 8 ); k 29 = 1.2727 x 10 6 
(2.2288 X 10®) ; and k 2 g = k 27 = k 39 = k 3 j = 0. 
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TABLE HI.- SUMMARY OF LONGITUDINAL FREQUENCIES AND DAMPING OF 
1/10-SCALE APOLLO/SATURN V MODEL 


S-IC, 

percent 

full 

Rigid body 
frequency, cps 

Mode 1 

Mode 2 



Mode 3 

Mode 4 

Mode 5 

f e , cps 

2? 

f 0 , cps 

f e , cps 

2? 


E 0 , cps 

f e , cps 

2? 

f 0 , cps 

f e , cps 

2? 

f 0 , cps 

f e , cps 

2? 


100 

3.5 

30.9 

— 

38.0 

40.2 

0.062 


40.3 

56.4 

0.025 

55.3 

71.3 

— 

70.8 

136.2 

0.030 

— 

75 

3.8 

41.2 

0.054 

— 

55.2 

.036 

j 

f 55.0 

— 

— 

— 

— 

— 

— 

— 

— 

-- 

50 

4.2 

42.0 

.068 

43.1 

57.1 

.026 

1 

1^58.2 

71.2 

.016 

70.8 

— 

-- 

— 

136.2 

— 

-- 

25 

4.8 

45.8 

.045 

— 

59.0 

.030 

1 

(”70.8 

86.6 

— 

— 

94.2 


— 

139.5 


— 

0 

5.7 

51.1 

.023 

53.3 

69.6 

.027 


1^74.7 

86.0 

.048 

78.0 

— 

— 

— 

135.1 

.055 



.. . 
























Figures.- 1/10- scale model. S- 1 1 stage fully assembled. 
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Figure 6.- I/10-scale model. S-ll simulated engines. 
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Figure 7.- 


I/10-scale model. S-IVB stage final assembly. 
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Figure 8.- 1/10-scale model. Simulated LM mounted in adapter structure. 
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Figure 9.- 1/10-scale model. Support cradle and cross beam 
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Figure 11.- 1/10-scale model. Typical accelerometer attachment. 
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Model station 


m 




in. 

m 

4l6.0 

10.566 

377.1 

9.578 

360.3 

9.151 

342.0 

8.687 

332.6 

8.448 

327.4 

8.316 

321.1 


312.2 

7.930 

275.6 

7.000 

266.0 

6.756 

254.5 

6.464 

239.5 

6.083 

183.4 

4.658 

167.1 

4.244 

15.6.4 

3.972 

146.9 

3.733 

84.7 

2.151 

64.2 

1.631 

33.0 

0.838 

20.5 

0.521 

12.4 

0.315 

—2.0 

-0.051 


Orientation 

Longitudinal 

Longitudinal 

Longitudinal 

LM Tri-axial 
Longitudinal 
Longitudinal 

Longitudinal 


Longitudinal 

Longitudinal 

Pitch and yaw — engine center line 
Longitudinal 


Longitudinal 

igitudinal-thrust structure center 
Pitch and yaw — engine center line 
Longitudinal 


— 84.7 2.151 Longitudinal 

— 64.2 I. 65 I Longitudinal 

33*0 O .838 Longitudinal 

— 20*5 0.521 Longitudinal-thrust structure center l: 

12.4 0.315 Longitudinal 

— —2.0 -0.051 Longitudinal — center line all engines 
Figure 12.- Schematic of transducer locations and orientations. 
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(e) S - 1C general layout 
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Figure 13.- Schematic of strain and pressure transducer locations and orientations in model S-IC stage. 
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Figure 15.- Lumped-parameter analytical model. 






















Normalized model station, x/L 
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Figure 27.- Experimental and analytical deflection shapes for second resonance of 100-percent full condition, 
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Figure 28.- Experimental and analytical deflection shapes for third resonance of 100-percent full condition. 
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Normalized model station, x/L 



Figure 29.- Experimental and analytical deflection shapes for fourth resonance of 100-percent full condition. 
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Figure 35.- Experimental and analytical deflection shapes for first resonance of 50-percent full condition. 
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Figure 38,- Experimental and analytical deflection shapes for third experimental resonance of 50-percent full condition. 
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Figure 43.- Experimental deflection shape for third resonance of 25-percent full condition. 
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Figure 50.- Experimental and analytical deflection shapes for third resonance of burnout condition. 
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Figure 54.- Static deflection shape for 400 lb (1780 N) longitudinal load at y- = 0.88. 
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Figure 57.- Repeating elements of S-IC LOX tank. 
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Figure 58.- Configuration and moment of inertia of one S-ll cross beam. 
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